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Abstract. Recent far-ultraviolet (FUV) absorption line measurements of low-redshift quasars have unveiled a 
population of intervening broad Hi Lya absorbers (BLAs) with large Doppler parameters {b > 40 kms~^). If 
the large width of these lines is dominated by thermal line broadening, the BLAs may trace highly-ionized gas 
in the warm-hot intergalactic medium (WHIM) in the temperature range T « 10^ - 10*^ K, a gas phase that 
is expected to contain a large fraction of the baryons at low redshift. In this paper we use a hydrodynamical 
simulation to study frequency, distribution, physical conditions, and baryon content of the BLAs at z « 0. From 
our simulated spectra we derive a number of BLAs per unit redshift of (dA/" /dz)BLA ~ 38 for Hi absorbers with 
log (A/'(cm"^)/6(kms"^)) > 10.7, fo > 40 kms"\ and total hydrogen column densities A''(Hll) < 10^°'^ cm~^. 
The baryon content of these systems is r2i,(BLA) = 0.0121 ftes which represents ~ 25 percent of the total 
baryon budget in our simulation. Our results thus support the idea that BLAs represent a significant baryon 
reservoir at low redshift. BLAs predominantly trace shock-heated collisionally ionized WHIM gas at temperatures 
log T ~ 4.4 — 6.2. About 27 percent of the BLAs in our simulation originate in the photoionized Lya forest (log 
T < 4.3) and their large line widths are determined by non-thermal broadening effects such as unresolved velocity 
structure and macroscopic turbulence. Our simulation implies that for a large-enough sample of BLAs in FUV 
spectra it is possible to obtain a reasonable approximation of the baryon content of these systems solely from the 
measured H I column densities and b values. 

Key words, methods: numerical - cosmology: miscellaneous - large-scale structure of Universe 



• I— I 1. Introduction of bright extragalactic background sources. Space-based 

observatories such as the Far Ultraviolet Spectroscopic 

H , Cosmological simulations predict that a large fraction Explorer (FUSE), the Hubble Space Telescope with its 

. , (~ 30 - 40 percent) of the baryons m the local Universe Space- Telescope Imaging Spectrograph (STIS), and the 

exist m the form of hot (T = 10^ - 10^ K), highly ion- chandra and XMM-Newton X-ray observatories have been 

ized intergalactic gas (e.g., Cen & Ostriker 1999; Dave used to the study frequency and distribution of highly ion- 

et al.2001), commonly referred to as the warm-hot m- .^^^ ^^^^^^ ^q^^^ q^j^^ q^j^) ^jjIM toward a 

tergalactic medium (WHIM). The WHIM is believed to ^^^^^er of low-redshift background sources (e.g., Tripp, 

emerge from intergalactic gas that is shock-heated to high g^^^g^ ^ j^^l^i^^ 2000; Savage et al. 2002; Fang et al. 2002; 

temperatures as the medium is coUapsmg under the action g^^^^ ^ Danforth 2005; Nicastro et al. 2005) . FUV obser- 

of gravity (Valageas, Schaeffer, & Silk 2002). With current ^^^^^^^ .^p^^ ^^^^ ^^^^^^ ^^^^^^^ q absorbers 

instrumentation, the most promising method to study the ^^^^^^ beheved to trace the WHIM at temperatures 

distribution and baryon content of the WHIM IS to search ^^^^^^ T = 3 x 10^ K) is ^(^(Ovi) > 0.0022 /lyo"' 

for absorption hues from neutral hydrogen and highly- (Danforth & ShuU 2005; Sembach et al.2004), but this 

ionized heavy elements m far-ultraviolet and X-ray spectra estimate depends critically on the assumed oxygen abun- 

dance, which is not well constrained. 

Send offprint requests to: P. Richter 

e-mail: prichter@astro.uni-bonn.de Apart from highly ionized oxygen, recent STIS ob- 

* Chandra Fellow servations by Richter et al. (2004, 2006) and Sembach et 

** Founded by merging of the Institut fiir Astrophysik al. (2004) suggest that WHIM filaments can be traced by 

und Extraterrestrische Forschung, the Sternwarte, and the thermally broadened H I Ly a absorption (see also Williger 

Radioastronomisches Institut der Universitat Bonn. et al. 2006). Although the vast majority of the hydrogen 
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Fig. 1. Three examples for BLA absorbers in our simulated spectra are shown. The panels show the logarithmic total 
hydrogen volume density, gas temperature, neutral hydrogen volume density, and normalized intensity for Hi Lya 
and Ovi A1031.9 absorption as a function of the radial velocity along each sightline. Densities n and temperatures T 
are in units cm~"^ and K, respectively. 



in the WHIM is ionized due to particle collisions and the 

ambient UV background radiation, a very small fraction 
(< 10~^, typically) of the hydrogen is expected to be neu- 
tral if the gas is in ionization equilibrium. For WHIM fil- 
aments with total gas column densities > 10-*^^ cm~^ the 
neutral hydrogen should give rise to weak Lya absorp- 
tion with H I column densities on the order of log N !v 13. 
Due to the high temperature of the gas (T = 10^ — 10^ 
K) and the effect of thermal line broadening, the Ly a ab- 
sorption is expected to be very broad and shallow with 
large Doppler parameters {b > 40 kms"""^). About fifty of 
these broad Lya absorbers (BLAs) have been identified 
in STIS low-z QSO spectra so far (Richter et al. 2006), in- 
dicating that these systems represent a significant baryon 
reservoir at low redshift. The STIS measurements imply 
a number of BLAs per unit redshift of (dAf /dz)Bi,A > 22 
for absorbers with log (iV(cm~^)/5(kms~-'^)) > 11.3 (de- 
tection limit for BLAs in STIS data with a signal-to- 
noise ratio of ~ 18 per 10 kms~^ wide pixel, hereafter 
referred to as (S/N)io), and a baryon mass density of 
176(BLA) > 0.0027 /i7o"^ The advantage of using the 
BLAs for an estimate of fib is that no assumption for the 
metallicity of the WHIM gas has to be made and there- 
fore also low- metallicity environments can be detected. 
Large uncertainties for f26(BLA) arise, however, from the 
unknown contribution of non-thermal line broadening ef- 
fects and the uncertain ionization fraction of the gas (see 
Richter et al.2006). 

In this paper we use artificial spectra generated from 
a cosmological simulation to study the frequency and 
baryon content of BLAs at low redshift. The aim of 
this paper is to test whether these simulations reproduce 
the observed spectral signatures of WHIM filaments (i.e., 
broad and shallow Hi Lya absorption) and to investi- 



gate the physical conditions in these systems. This paper 
is organized as follows: in Sect. 2 we explain the simula- 
tion method; in Sect. 3 we describe the spectral analysis 
of the artificial spectra; physical properties of the BLAs 
arc discussed in Sect. 4; in Sect. 5 we estimate the baryon 
content of the BLAs in our simulation and compare the 
results with recent observations; in Sect. 6 we compare the 
relationship between BLAs and O VI absorbers; finally, we 
summarize our study in Sect. 7. 

2. Simulation metliod 

For our study we use the output of a numerical simulation 
that was part of an investigation of intervening O vi ab- 
sorption arising in WHIM filaments (Fang & Bryan 2001), 
and which is based on a grid-based adaptive mesh refine- 
ment method (Bryan 1996; Norman & Bryan 1999). The 
generated cube is 20 Mpc on a side and has a gas and 
dark matter mass resolution of 6 x 10^ and 5 x 10^ Mq, 
respectively. The maximum resolution is 9.8 kpc (smallest 
grid cell) . The initial density field was drawn from an adi- 
abatic cold dark matter power spectrum as approximated 
by the relation given in Eisenstein & Hu (1998). The sim- 
ulation is based on a cosmological model with = 0.4, 
fib = 0.0473, Oa = 0.6, and h = 0.65 (in units of 100 
kms""'^ Mpc""'^). The inital redshift is z = 50 and the sim- 
ulation stops at z = 0.2. Radiative heating and cooling 
of the gas is not included in the simulation, so that the 
gas temperature in the unshocked, low-density regions is 
significantly underestimated. For these regions we there- 
fore approximate the temperature of the gas using the 
temperature-density relation T = Tg (1 + S)''^^ found for 
the Lya forest, where Tq « 5000 K, 7 « 1.4 (e.g., Zhang 
et al. 1997) and 6 is the overdensity of the gas. Note that 
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Fig. 2. Cumulative number of BLAs per unit redshift, log 
{dJV/dz)BLA, as a function of the absorption strength, 
log (iV(cm~^)/&(kms~^)). The total sample is shown in 
gray, the high-quality sample is plotted in black. The sin- 
gle data point indicates the value for {dJV/dz)BLA from 
FUV absorption line measurements with STIS (Richter et 
al. 2006). 



previous studies indicate that radiative heating/ cooling is 
not expected to significantly influence the evolution of the 
WHIM gas at low redshifts (see, e.g., Dave et al. 2001). 
However, we are planning to investigate this issue in more 
detail in a future paper. A density-dependent metallic- 
ity model was adopted, based on the results from Cen 
& Ostriker (1999). More detailed information about the 
simulations can be found in Fang & Bryan (2001) and 
Norman & Bryan (1999) and references therein. The ion- 
ization state of the gas in the simulation has been mod- 
elled using the code CLOUDY (version 90.04; Ferland et 
al. 1998) , which includes both coUisional ionization and 
photoionization. A grid of temperatures (log T = 3 — 7) 
and densities (log nn = —8 to —2) is calculated from 
CLOUDY assuming an ionizing background spectrum at 
z = 0.2 based on the model of Haardt & Madau (1996) 
together with a mean specific intensity at the Lyman limit 
of = 2 X 10-23 ergs-^Hz-^sr-i. 

For the analysis of the simulated region we have drawn 
3000 random lines of sight (LOS) through the simulated 
box and have created artificial spectra for Hi and Ovi 
following the method described in Zhang et al. (1997). The 
total redshift path of these 3000 sightlines is Az = 20.0. 
The spectra then have been fitted using the automatted 
line fitting algorithm AutoVP (Dave et al. 1997), which 
delivers column densities, N, and Doppler parameters, b, 
for the fitted spectral features. 



3. Hi Spectral analysis 

3.1. BLA sample selection 

The fitting of the 3000 LOS with AutoVP results in the 
detection of 31,340 Lya absorbers. The vast majority of 
these lines belong to the photoionized Lya forest. 2671 of 
these lines have 6 > 40 kms"^ and thus represent BLA 
candidates that may trace collisionally ionized gas at tem- 
peratures T > 10^ K arising in WHIM filaments. We have 
applied a number of selection criteria to identify the best 
BLAs candidates in our sample. First, we have excluded 
broad features that are part of complex multi-component 
absorbers, as the determination of line centers, b values, 
and column densities for these systems is afilicted with 
large systematic uncertainties. To exclude these crowded 
regions, we take into account only those BLA candidate 
lines whose nearest neighbouring absorption line is at least 
2a away, where a is the standard deviation of a Gaussian- 
shaped line defined as: 



a = 



FWHM 
2.354 



= 0.705 6. 



(1) 



In addition, wq have excluded very weak BLA lines by 
constraining our BLA sample to H i absorbers that have an 
absorption strength of log {N{cm~^)/b{kms~^ j) > 10.7. 
This limit corresponds to a minimum absorption depth of 
~ 10 percent. Note that lines with log (N/b) < 10.7 do not 
significantly contribute to the baryon content of the BLAs; 
moreover, they would clearly fall below the detection limit 
in FUV spectra that have signal-to- noise ratios (S/N)io < 
50. After applying the above selection criteria we are left 
with a sample of 770 BLA candidate lines. 

BLAs arising in the WHIM probably do not repre- 
sent single-component, isothermal gaseous structures, but 
rather are expected to be multi-phase systems with den- 
sity substructure and internally varying ionization con- 
ditions. The method of using the BLAs to measure the 
baryon content of the WHI^I is based on the idea that the 
shock-heated regions contain tlie dominant fraction of the 
baryons in these systems, and that these regions also de- 
termine the Ly a line width due to their high temperature 
and the resulting thermal line broadening. In an ideal case, 
one thus would expect to observe a single-component, 
Gaussian-shaped broad absorption feature with a line- 
width that serves directly as a measure of the tempera- 
ture of the gas (see also Sect. 3.4). However, if substan- 
tial substructure is present in a BLA, the thermal line- 
broadening is competing with the broadening of the ab- 
sorption due to unresolved radial-velocity components. To 
distinguish between thermal broadening and non-thermal 
broadening it is therefore important to look for evidence 
of velocity-component structure in the BLA candidate 
lines (see also Richter et al. 2006). We have done a vi- 
sual inspection of our 770 BLA candidates to evaluate the 
quality of the line fitting (as done by AutoVP) and to 
look for evidence for sub-component structure. 321 BLAs 
in our sample are almost perfectly Gaussian shaped and 
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thus represent excellent candidates for thermally broad- 
ened WHIM absorbers. The remaining 449 BLAs show 
evidence for sub-component structure that has been ig- 
nored by the automated line fitting procedure. It could 
be that a large fraction of these broad systems may not 
be related to the WHIM but rather represent photoion- 
ized multi-component Lya forest lines. Note that in real 
FUV spectra with S/N ratios (S/N)io < 50 these multi- 
component systems would not be distinguishable from the 
single-component cases. To compare our results with ob- 
servational data it is therefore important to assess the 
properties of all BLA candidate lines in our sample. For 
the following analysis we distinguish between the high- 
quality BLA sample (321 absorbers) and the total BLA 
sample (770 absorbers). 

Examples for BLAs along different lines-of-sight (LOS) 
in our simulated spectra are presented in Fig. 1, where we 
have plotted the total hydrogen volume density, the gas 
temperature, the neutral hydrogen volume density, the 
normalized Hi Lya A1215.7 intensity, and the normal- 
ized intensity of the Ovi A1031.9 line as a fmiction of the 
restframe velocity for three BLAs from the high-quality 
sample. As clearly visible, the Hi Lya absorption in these 
systems is broad and symmetric. The left example (LOS 
2410) shows a WHIM BLA with 6 « 59 kms'^ and log 
iV(H i) = 13.37; the density and temperature distribution 
is very simple in this system. The width of this BLA line 
is governed mainly by thermal broadening of the WHIM 
gas that has log T « 5.3. The middle panel (LOS 1166) 
shows a thermally-broadened BLA with b w 129 kms~^ 
and log iV(Hi) = 13.28. Also this system has a simple 
single-component structure, but the WHIM temperature 
is with log T ss 5.9 somewhat higher than in the left ex- 
ample, thus resulting in an Hi Lya line that is broader 
(and shallower). The right panel (LOS 276) shows an ex- 
tremely broad BLA system with b « 165 kms""'^ and log 
Af(Hi) = 13.70. Here, the BLA system is characterized by 
a complex density distribution that spans several hundred 
kms^^ (see upper panel). The BLA has a high tempera- 
ture of log T ft: 6.0, but in this case, the extremely large 
width of the line is not dominated by thermal broadening 
but by the c;oniplex vclocdty structure. The shape of the H I 
Ly a absorption slightly deviates from a perfect Gaussian, 
indicating the presence of velocity sub-components in this 
line. Such small asymmetries would remain unnoticed in 
real spectra where the S/N is limited. 



3.2. BLA number density 

Fig. 2 shows the cumulative number of BLAs per unit red- 
shift, log [dM/dz), plotted against the BLA absorption 
strength, log {N/b). The total sample (770 systems) is 
shown in gray, the high-quality sample (321) is plotted 
in black. For 6 > 40 kms~^ and log (N/b) > 10.7 our 
simulation predicts that {dJV / dz)B'LA ~ 39 (total sam- 
ple) and (cLV /dz)BLA ~ 16 (high-quality sample). For log 
{N/b) > 11.3 the BLA number density in the total sam- 
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Fig. 3. Hi column density and b value distribution result- 
ing from the Voigt fit of the BLA lines in the simulation. 
The gray line shows the distribution of the total sample 
(770 systems), the black line shows the distribution for 
the restricted sample (321 systems). Column densities are 
in units cm~^, b values are in units kms~^. 



pie reduces to {dM /dz)BLA ~ 19. This is very similar to 
the value obtained by recent STIS observations of BLAs 
towards low-redshift quasars {{dj\f / dz)BLA = 22 for reli- 
ably detected systems; Richter et al. 2006). The value for 
{dM /dz)BLA from the STIS measurements is indicated in 
Fig. 2. 

3.3. Hi column densities 

Fig. 3 shows the distribution of logarithmic H i column 
densities and logarithmic H i 6 values, as obtained from 
the profile fitting of the BLA candidate lines. The total 
sample is shown in gray, the high-quality sample is shown 
in black. The Hi column density distribution (Fig. 3, up- 
per panel) shows that the majority of the BLA systems 
(~ 92 percent) have relatively low Hi column densities 
with log 7V(Hi) < 14 and a peak near log 7V(Hi) « 13. 
Note that for column densities lower than that our sam- 
ple is incomplete due to the exclusion of systems with log 
{N/b) < 10.7 (see Sect. 3.1). The true Hi column density 
distribution of the BLAs most likely further rises towards 
lower column densities. The distribution of BLA H i col- 
umn densities in our simulation is qualitatively very simi- 
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lar to the one observed in STIS spectra of low— z quasars 
(Richter et al. 2006; their Fig. 5). 

3.4. Hi b values 

The BLA b value distribution is plotted in the lower panel 
of Fig. 3. Some BLAs in our simulation have b values larger 
than 200 kms~^, but the majority of the systems have rel- 
atively moderate line widths with b values ranging from 
40 to 100 kms-i (log 6 = 1.6 - 2.0). The 6 values of the 
BLAs are assumed to be composed of a thermal compo- 
nent, 6th, and a non-thermal component, fenon-th, in the 
way that 



^non — 



th 



(2) 



The non-thermal component may include pro- 
cesses like macroscopic turbulence, unresolved velocity- 
components, and others (see Richter et al. 2006 for a de- 
tailed discussion). The contribution of the thermal com- 
ponent to b depends on the gas temperature: 



bth = 



2kT 



0.129i 



■ kms 



(3) 



where k is the Boltzmann constant, m is the particle 
mass, and A is the atomic weight. For the shock-heated 
WHIM gas with log T > 5 one thus expects &th > 40 
kms""*^. The non-thermal broadening mechanisms are ex- 
pected to contribute to some degree to the total b values 
in WHIM absorbers (see Richter et al. 2006), so that the 
measured b value of a BLA provides only an upper limit for 
the temperature of the gas. If we compare the b value dis- 
tribution of the total BLA sample in our simulation (lower 
panel in Fig. 2, gray line) with that of the high-quality 
sample (black line), we see that the fraction of absorbers 
with log b < 1.8 {b < 63 kms^^) is very large in the total 
sample (^ 70 percent), while for the high-quality sample 
it is somewhat smaller (^ 59 percent). This implies that 
for low b values the contamination from multi-component 
Lya absorbers (i.e., systems for which 6non-th is dom- 
inant) in the total BLA sample is substantial (see also 
Sects. 3.1, 4.1, and 5.1). 

4. Physical properties of BLAs 

4.1. Temperature distribution 

We have determined for each BLA in our sample a charac- 
teristic temperature by calculating the mean of the tem- 
peratures in the cells that contribute to the Hi Lya ab- 
sorption, weighted by Hi optical depth in each cell. The 
Hi weighting is a crucial point, since it is the tempera- 
ture in the cells with the largest Hi optical depths that 
determines the thermal broadening (and thus the width) 
of the Ly a absorption. In Fig. 4, upper panel, wc show the 
distribution of log T for the total sample (gray) and the 
high-quality sample (black) . The BLA temperature distri- 
bution is bimodal, with most BLAs spanning a range of 




17 18 19 20 
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Fig. 4. Distribution of temperatures (T in K), hydrogen 
volume densities [n in cm"'^), and total hydrogen column 
densities {N in cm^^) for the total BLA sample (770 ab- 
sorbers; shown in gray) and the high-quality sample (321 
absorbers; shown in black). 



4.4 < logT < 6.2 with a peak near log T « 5.2, thus typ- 
ical for what is expected for the low-temperature WHIM 
(e.g., Dave et al. 2001). A second sharp peak is visible 
at lower temperatures near log T = 3.8, at temperature 
that is characteristic for the photoionized Lya forest (e.g., 
Dave & Tripp 2001). At such low temperatures, thermal 
broadening contributes only with 6th ~ 10 kms~^ to the 
total line width of these systems (see equation (3)). The 
measured large width of these lines (6 > 40 kms~^) there- 
fore must be determined by non-thermal processes, such 
as unresolved component structure or macroscopic turbu- 
lence. 

The BLA temperature distribution clearly separates 
low-tcmpcrature, photoionized H I absorbers with log T < 
4.3 (i.e., broad Lya forest lines) from shock-heated, high- 
temperature WHIM absorbers that have log T > 4.3 (in- 
dicated by the dotted line in Fig. 4, upper panel). While 
for the high-quality sample only about eight percent of 
the BLAs have temperatures log T < 4.3, the fraction of 
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BLAs with log T < 4.3 in the total sample is ^ 27 per- 
cent. Nearly all (198 out of 209) of these low-temperature 
BLAs in the total sample have b values < 63 kms^^ (log 
b < 1.8). This shows that BLA absorbers with b values 
< 63 kms~^ in the total BLA sample are substantially 
contaminated by photoionizcd Lya forest lines (see also 
Fig. 3 and Sect. 3.1), while for BLAs with b > 63 kms~^ 
the collisionally ionized WHIM absorbers dominate the 
BLA population. 

Note that due to the relatively small box size and the 
missing large-scale power in our simulation the charac- 
teristic temperatures of WHIM filaments and BLAs may 
be slightly underestimated. We are planning to examine 
this issue in more detail in a future paper using a higher- 
resolution simulation. 

4.2. Density distribution 

The distribution of logarithmic hydrogen volume densi- 
ties (riH in units cm~^) in our two BLA samples is shown 
in Fig. 4, middle panel. In a similar way as for the tem- 
peratures, we have calculated for each BLA nn from the 
mean of the hydrogen volume densities in the cells that 
contribute to the H i absorption, weighted by the H i opti- 
cal depth in each cell. The distribution of log tih spans a 
large range between —3.6 and —7.2. Similar to the temper- 
ature distribution, the distribution of log nn is bimodal 
with one peak near log nu ~ — 5.4 and a second one near 
log Till « —6.7. Again, this bimodal distribution sepa- 
rates the low-density, photoionized Lya forest lines (log 
riH < —6.2, typically) from the higher-density WHIM ab- 
sorbers (log TiH > —6.2, typically). The above density 
peaks near log tih = —6.7 and —5.4 correspond to over- 
densities of ph /ph ^ 1 for the low-density Ly a forest and 
Ph/ph ~ 20 for the WHIM absorbers, thus similar to what 
has been found in other simulations (e.g., Dave & Tripp 
2001; Dave et aL2001). 

4.3. Ionized hydrogen column densities 

The column density of ionized hydrogen in each BLA 
(A^(Hii) « -^^(H)) can be estimated directly from our sim- 
ulation, i.e., by integrating the hydrogen volume density, 
riH, over the cells that are located within the 2a velocity 
range of each BLA (see equation (1)). In the lower panel 
of Fig. 4 we show the distribution of log A^(H ii) for the 
total sample (gray) and the high-quality sample (black). 
Most BLAs have ionized (=total) hydrogen column den- 
sities in the range log A^(H ii) = 18 — 20 with a peak near 
log A''(Hii) = 18.8. Less than seven percent of the BLAs 
have Hii column densites larger than 10^° cm~^, but note 
that these systems contain a substantial fraction of the to- 
tal gas column density in the BLAs (^ 57 percent in the 
total sample, ~ 51 percent in the high-quality sample). 87 
percent of the broad Lya forest lines with log T < 4.3 
have Hii column densities log A''(Hii) < 18.5. Broad ab- 
sorbers that arise in the WHIM therefore dominate the 
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Fig. 5. Upper panel: distribution of ionization fractions in 
the total BLA sample (gray) and the high-quality sam- 
ple (black). Lower panel: the hydrogen ionization fraction 
for each cell in our simulation (photoionization-|-collisional 
ionization), log /h = log(^^H/^^HI), is plotted as a function 
of the gas temperature, log (r(K)). The light gray shaded 
area indicates cells in the density range log nn = —5 to 
—7. The dark gray shaded area refers to cells that have log 
riH = —5.3 to —5.6, thus a density range that is character- 
istic for BLAs. The solid black line shows the best linear 
fit of the data in the BLA density range, the dashed black 
curve shows a function that approximates log f{T) for a 
pure coUisional ionization model. 

distribution of the BLA H ii column densities in the range 
log 7V(Hii) = 18.5-21.0. 

4.4. Ionization fraction 

In the upper panel of Fig. 5 we have plotted the dis- 
tribution of the ionized hydrogen ionization fractions, 
iV(Hll)/Af(Hi), for our two BLA samples. The logarith- 
mic ionization fractions in both samples typically range 
between 5 and 7 with a peak near 5.6. 

In the following we want to explore in more detail 
the individual ionization mechanisms (collisional ioniza- 
tion and photoionization) and their role in the observed 
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A'^(Hii)/A^(Hi) distribution. For pure coUisional ioniza- 
tion and in ionization equilibrium, the hydrogen ioniza- 
tion fraction, log /h = log (nH/jT-Hi); is governed solely 
by the gas temperature (e.g., Sutherland & Dopita 1993). 
The ionization fraction then can be approximated via 

log/H(r) « -13.9 + 5.41ogT-0.33(logT)2. (4) 

However, as the densities in the WHIM absorbers are 
low (log riH < —5, typically; see Sect. 4.2), it is likely that 
photoionization contributes to the ionization fraction of 
the gas. To assess the temperature dependence of /h in 
our combined coUisional+photoionization model we have 
plotted in Fig. 5, lower panel, for each of the ^ 1.5 x 10^ 
cells along our 3000 LOS the ionization fraction, log /h, 
as a function of the gas temperature, log T. The light gray 
shaded area shows cells with hydrogen volume densities in 
the range log uh = —5 to —7. The dark gray shaded area 
shows cells with log nu = —5.3 to —5.9, thus in a density 
range that is characteristic for the shock-heated BLAs (see 
Fig. 3). The dashed black curve indicates the Jh—T rela- 
tion for pure collisional ionzation (equation (4)). As clearly 
visible, all points lie above the relation for coUisionally ion- 
ized gas, demonstrating that photoionization is important. 
For the density range that is typical for the BLAs in our 
simulation (log nn = —5.3 to —5.9; dark gray shaded area) 
the combined collisional-|-photoionization model suggests 
that log /h is roughly linearly related to log T, so that 

log/H(r)«-0.75 + L251ogT. (5) 

This approximation is shown in Fig. 5, lower panel, as 
a solid black line. 

4.5. Non-thermal line broadening in BLAs 

In Fig. 6, upper panel, we have plotted the logarithmic b 
values from the line fitting of our BLA sample against log 
6th, which represents the thermal contribution to b for each 
BLA, as calculated from the BLA temperatures (Sect. 4.1) 
and the inversion of equation (3). The gray dots show the 
data points for our total BLA sample, the black dots show 
the data points for the; high-quality sample. For pure ther- 
mal broadening, b = bth (equation (2)). This relation is 
shown as solid black line in the upper panel of Fig. 5. As 
expected from the BLA temperature distribution (Fig. 4), 
there are a number of low-temperature, photoionized sys- 
tems that have 6th ^ b because of non-thermal line broad- 
ening effects; these systems cluster near log 6th = 1- The 
BLAs that refer to the shock-heated WHIM roughly fol- 
low the slope of the relation 6 = 6th, but are systematically 
offset by a few percent. A statistical comparison between 
6 and 6th for log 6th > 1-5 implies that for both the total 
sample and the high-quality sample the relation 6=1.1 6th 
gives the best fit to the data. This means that non-thermal 
broadening mechanisms contribute (on average) with ~ 10 
percent to the total 6 value of a WHIM BLA in our sim- 
ulation. From equation (2) then follows that (on average) 
^non-th ~ 0.5 6th for thcsc systems. Note that there are a 



7 




17 18 19 20 21 22 



log N(Hll)3i^ 

Fig. 6. Upper panel: The total 6 values (6 in units km s~^) 
for the BLAs as obtained from the profile fitting is plotted 
against the thermal component, 6th, as calculated for each 
BLA from the Hi optical depth weighted temperature 
in each system (Sect. 4.1). The solid black line indicates 
^tot = ^th- Lower panel: comparison between the total hy- 
drogen column densities {N in units cm~^) estimated in- 
directly from the measured line widths and H I column 
densities (log iV(Hii)fit; see Sect. 4.6) and derived directly 
from the simulation (log A/'(Hii)sii„; Sect. 4.3). The black 
solid line indicates log A''(Hii)fit =log A''(Hii)sim. 



number of systems that have 6 < 6th- Although this ap- 
pears nonsensical at first, an inspection of these systems 
suggests that their absorption profile is dominated by a 
cooler gas component that has a relatively large H i optical 
depth. The fitted 6 value in these systems only reflects the 
temperature of this cooler gas component, but is not good 
measure for the temperature of the hot gas, whose absorp- 
tion is hidden in a broad wing that is ignored by Auto VP. 
Although 6th is based on the optical-depth weigthed tem- 
perature of the gas in our BLAs, 6 < 6th can occur in a 
few systems where this effect is particularly pronounced. 
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4.6. Indirect estimate of the ionized gas content in 
BLAs 

The combination of the equations (2), (3), (5), and b = 
1.1 bth allows us to estimate the ionization fraction, /h, in 
a thermally broadened BLA indirectly from its measured 
b value. Together with the measured Hi column density 
one thus can estimate the total hydrogen column density 
in a BLA via Af(Hii) = /h-/V(Hi). Combining all these 
equations gives 

log TV (H II) « logiV(HI) + 1.25 log (0.816^) + 1.47. (6) 

Since b and iV(Hi) are the only parameters of BLAs 
that can be measured in real FUV spectra, this indirect 
estimate of Af(Hii) is the only possible method to ob- 
servationally constrain the total gas budget and thus the 
baryon content of BLAs in the low-redshift Universe (see 
Richter et al. 2006). With our simulation we now want 
to examine the reliability of this indirect method for our 
sample of BLA absorbers. In Fig. 6, lower panel, we have 
plotted A^(Hii)fit against 7V(Hii)sini for the total sam- 
ple (gray dots) and the high-quality sample (black dots). 
A'^(Hll)fit is the hydrogen column density derived indi- 
rectly via equation (6) and A^(Hii)sim is the hydrogen col- 
umn density obtained directly from the simulation using 
the cell integration (Sect. 4.3). The black solid line indi- 
cates the ideal case, iV(Hii)fit = A^(Hii)sim- Although the 
data points show a significant scatter around that line, the 
plot implies that (for a large enough sample) the indirect 
method indeed serves as a valuable tool to roughly esti- 
mate the total hydrogen column densities in BLAs. For 
the total sample the mean scatter is 0.25 dex, thus larger 
than for the restricted sample (0.18 dex) owing to the 
fact that non-thermal line broadening processes in these 
systems cause a considerable uncertainty for the indirect 
estimate of A/'(Hii). 

5. The baryon content of BLAs 

5.1. The baryon content of BLAs in the simulation 

From the cell integration (direct method; Sect. 4.3) we de- 
rive a total hydrogen column density of 2.31 x 10^^ cm~^ 
for the total BLA sample (770 systems) and 9.67 x 10^^ 
cm~^ for the high-quality sample (321 systems). The cos- 
mological mass density of the BLAs in terms of the current 
critical density Pc is given by 

0,(BLA) = ™^ 7V(HII)totai AX-\ (7) 

Pc C 

where = 1.3, mn = 1.673 x IQ-^^ kg, Hq = 65 
kms~^ Mpc~^, and Pc = 3Hq ^/8ttG. To compare our re- 
sults with recent observations we assume at this point 
that the comoving path length for our BLA sample is 
AX = Az = 20.0 (i.e., we assume that qo = 0). Using 
equation (7), the values for A^(Hii)totai listed above cor- 
respond to baryon densities of fi(,(BLA) = 0.0167 /i65~^ 
and 0.0070 /i65 respectively. From this follows that the 



BLAs trace a considerable fraction of the baryons in our 
simulation (~ 35 percent in the total sample, ~ 15 percent 
in the high-quality sample). Note that the contribution to 
r2f,(BLA) from photoionized absorption systems with log 
T < 4.3 (Lya forest) is negligible (~ 2 percent in the to- 
tal sample and < 1 percent in the high-quality sample) . In 
the upper left panel of Fig. 7 we show the distribution of 
baryons as a function of the total hydrogen column den- 
sity for the direct method. The gray line refers to the total 
sample, the black line indicates the high-quality sample. 
Values for log A/'(Hii) are given in 0.5 dex- wide bins and 
the values for f75(BLA) indicate the baryon content per 
bin. The lower left panel of Fig. 7 shows the number of 
BLAs in each column-density bin. From this comparison 
we learn that BLA systems with log A''(H ii) = 19.5 — 20.5 
contribute most to fib (BLA) (~ 56 percent), while they 
contribute with only ~ 16 percent to the total number of 
BLAs. 

5.2. The baryon content of BLAs derived from the 
indirect estimate 

Using the indirect method for the determination of 
7V(Hii) (equation (6)), we obtain values of r2b(BLA) = 
0.0465/165"^ (total sample) and 0.0070 /les"^ (high- 
quality sample). Thus, while direct and indirect method 
give identical results for the high-quality sample, the in- 
direct method vastly overestimates f2f,(BLA) for the total 
sample. In the right panel of Fig. 7 we plot (in a way sim- 
ilar as for the direct method) the distribution of rif,(BLA) 
and BLA numbers for bins of log A'^(Hll). When com- 
pared to the left panel, it immediately becomes clear that 
the indirect method wrongly suggests the presence of a 
few (three) BLA systems with very large total hydrogen 
column densities (log iV(Hii) > 21.5) that would contain 
more baryonic matter than all the other BLAs together. 
This is also seen in the lower panel of Fig. 6, where most of 
the few points with log A^(IIll)sii„ > 20.5 clearly lie above 
the expected relation iV(IIii)iit = A^(II ii)sii„. A visual 
inspection of these BLA systems shows that all but one 
systems are multi-component absorbers whose widths are 
determined mostly by non-thermal line broadening. The 
velocity-component structure is not properly accounted 
for in the fit by AutoVP, resulting in large discrepancies 
between iV(IIii)fit and A^(IIii)si,„ for these systems. As 
these high-column density absorbers dominate f2b(BLA), 
the overestimate of A^(Hii) for these systems with the in- 
direct method together with the low-number statistics in 
the high-column density regime introduces a substantial 
uncertainty for the (indirect) determination of f2b(BLA). 
The discrepancy in Jib (BLA) between the direct and in- 
direct method becomes less severe if we ignore these few 
high-column density systems. If we restrict our total BLA 
sample to values of log iV(Hii)fit < 20.5 (754 systems), 
we derive values of f7;,(BLA) = 0.0121 /igs ~^ (cell integra- 
tion) and 0.0136 /i65~^ (indirect method), thus a much 
better agreement between these two methods. In Table 1 
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Fig. 7. Distribution of f2(,(BLA) and BLA numbers for bins of log A^(Hii), where N is in units cm~^. The left panel 
shows the distribution based on the direct method to estimate Af(Hii) (cell integration), the right panel shows the 
distribution based on the indirect method (equation (6)). The bin size for log A''(Hii) is 0.5 dex. 



we compare values for fit(BLA) from the direct and the 
indirect method for a number of different selection crite- 
ria. From this follows that the agreement between direct 
and indirect method actually is quite good if one ignores 
the very broadest BLAs that have b values larger than 200 
kms~^. 

Following the statistics in our simulation, the best 
range for the analysis of BLA systems in future high- 
resolution, high S/N FUV spectra is 6 = 60 - 200 kms'^ 
For this range, our analysis yields {dJ\f / dz)-BhA ~ 13 and 
values for r2(,(BLA) of 0.0104 /igs (cell integration) and 
0.0111 /i65~^ (indirect method). Thus, while the agree- 
ment between direct and indirect method is excellent for 
these b values, the contamination of the BLA number den- 
sity by photoionized Lya forest lines with log T < 4.3 is 
negligible (~ 5 percent). 

5.3. Comparison with FUV observations 

Richter et al. (2006) have studied properties and baryon 
content of BLAs in a sample of four low-redshift quasars 
using STIS high-resolution spectra. Along a total red- 
shift path of Az = 0.928 they detect 49 BLA candidate 
systems above a sensitivity limit of log (N/b) « 11.3. 
However, only 20 of these candidates with b < 100 kms^^ 
are considered as secure detections. From these 20 BLAs, 
Richter et al. (2006) obtain values of (cW /dz)BLA = 22 
and rif,(BLA) = 0.0025 /igs It is important to note that 
their estimate of fib (BLA) is based on the assumption that 
the gas is in coUisional ionization equilibrium. Yet, our 



simulation shows that photoionization is important when 
it comes to the estimate of total hydrogen column densi- 
ties and f2b(BLA). If we take the 20 BLA candidates (se- 
cure detections) in the STIS data presented by Richter et 
al. (2006) and transform their measured Hi column densi- 
ties and 6 values (their Table 4) into total gas column den- 
sities using equation (6) (i.e., now taking both collisional 
ionization and photoionization into account), we obtain a 
total hydrogen column density of 3.45 x 10^" cm~^, corre- 
sponding to ili){BLA) = 0.0047 /i65 This estimate lies 
above the original estimate of 0&(BLA) = 0.0025 /igs ~^ by 
Richter et al., suggesting that the consideration of pho- 
toionization is important for the estimate of the baryon 
content of BLAs. If we restrict our total BLA sample in 
our simulation to log (N/b) « 11.3, log A^(Hii) < 20.5 
and b < 100 kms~^, we obtain {dJ\f/dz)BLA = 18 and 
12f,(BLA) = 0.0058 /i65 Thus, despite the small-number 
statistics for BLAs in the STIS data, observations and sim- 
ulations are in good agreement with each other. 



Our study suggests that the study of BLAs in low- 
redshift quasar spectra represents a useful method to 
study the distribution and baryon content of the WHIM. 
However, a number of BLA candidates significantly larger 
than currently available is desired to constrain the prop- 
erties of these absorbers on a statistically scicure basis. For 
this, future FUV measurements of BLAs at high spectral 
resolution and sufficient S/N will be required. 
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Table 1. Distribution of the baryons in BLAs with log (7V(cni~^)/6(kms^^)) > 10.7 



Selection 


N 


dU/dz 


f2j,(BLA)Bini /l65 


f26(BLA)fit /i65 


Total sample 


770 


39 


0.0167 


0.0465 


High-quality sample 


321 


16 


0.0070 


0.0070 


6 = 40 - 100 kms"' 


700 


35 


0.0076 


0.0102 


&= 100-200 kms"^ 


64 


3 


0.0064 


0.0068 


6 > 200 kms"^ 


6 


< 1 


0.0027 


0.0295 


b = 60 - 200 kms"^ 


268 


13 


0.0104 


0.0111 


log iV(Hll)flt < 20.5 


754 


38 


0.0121 


0.0136 


log Ar(Hii)fit < 20.0 


706 


35 


0.0076 


0.0077 


logiV(Hii)fit < 19.5 


597 


30 


0.0038 


0.0035 


log Ar(Hii)flt < 19.0 


409 


20 


0.0012 


0.0010 



6. The kinematic relation between BLAs and Ovi 
absorbers 

A considerable fraction of the BLAs in our sample show 
associated O vi absorption within a velocity radius of 50 
kms"'^. Number distribution and equivalent- width distri- 
bution of the O VI absorbers are discussed in detail in Fang 
& Bryan (2001). Here we focus on the kinematic relation- 
ship between broad Ly a absorption and associated O VI 
A1031.9 absorption in our simulation. Recent STIS obser- 
vations of BLAs and associated Ovi absorbers (Richter 
et al. 2006) show that many of the BLA/Ovi absorber 
pairs have velocity offsets between the Hi and the Ovi 
absorption. In addition, measured Ovi line widths often 
are inconsistent with those of the broad Hi absorbers (if 
a single-component absorption is assumed). These mea- 
surements imply that sub-component structure is present 
and/or that the broad Lya absorption and the Ovi ab- 
sorption do not trace the same gas phase. We now want 
to study these interesting aspects in our artificial spectra. 

From the 770 BLA systems in our sample, 274 systems 
have O vi absorption that is accurately fitted by AutoVP. 
Note that many of these Ovi absorbers are composed of 
several sub-components. For this analysis we consider only 
the dominating Ovi sub-component for each BLA with 
the highest Ovi optical depth. In Fig. 8, upper panel, we 
have plotted the distribution of the velocity offsets be- 
tween the line centroids of the H I and the O VI absorption, 
Av — |whi — woviI- O VI absorbers that refer to BLAs from 
the total sample (274 systems) are plotted in gray, Ovi 
systems that refer to BLAs in the high-quality sample (145 
systems) are plotted in black. For the total sample the me- 
dian offset between H l and O vi is 6.3 km s~ ^ , while for the 
high-quality sample the median offset is 5.3 kms~^. This 
implies that small velocity offsets between the H i Ly a ab- 
sorption and the O vi absorption are common. A possible 
explanation for such velocity offsets is that some of the 
gas filaments that give rise to broad H i and O vi absorp- 
tion have a relatively complex (non-Gaussian) density and 
temperature distribution (see, e.g., Fig. 1, right panel). 
While the BLAs trace the H i optical-depth weighted den- 
sity distribution of the filaments at all temperatures, O VI 
absorbers preferentially arise from higher- density regions 
in a small temperature range around 3 x 10^ K. Therefore, 



BLAs and Ovi absorbers most likely originate in some- 
what different gas phases and regions within a filament. 
This would naturally explain the observed small velocity 
offsets between these two species in both simulations and 
observations. 

In the lower panel of Fig. 8 we compare H i and O vi 
h values for the total sample (gray dots) and the high- 
quality sample (black dots). If Hi and O VI would trace the 
same gas phase of a purely thermally broadened WHIM 
absorber one would expect that &(Hi) = 4 6(Ovi) owing 
to the 16 times larger weight of oxygen (see equation (2)). 
Most of the BLA/Ovi absorber pairs lie above this ex- 
pected relation (in Fig. 8, lower panel, shown as solid line) 
and the data points show a large scatter for both the total 
and the high-quality sample. The plot indicates that the h 
values of H i and O Vi are only very poorly related to each 
other. This is not surprising in view of the fact that for 
most BLAs the absorption of H i and O vi is offset by a 
few kms~^ (see above). The Ovi absorbers have widths 
significantly larger than what is expected for pure ther- 
mal broadening in these systems, suggesting that most of 
the Ovi 6 values are determined predominantly by non- 
thermal processes. 

7. Summary 

In this paper we have used artificial spectra generated 
from a hydrodynamical simulation to study frequency, 
distribution, physical conditions, and baryon content of 
broad Lya absorbers (BLAs) in the low-redshift inter- 
galactic medium. Along a total redshift path of Az — 20 
we find 770 BLA systems with > 40 kms^^ above a 
detection threshold of log (A^(cm-2)/&(kms~^)) > 10.7, 
implying a number density of BLA systems per unit red- 
shift of {dAf / dz)BLA = 39 for this sensitivity limit. From 
a spectral analysis with automatted Voigt-profile fitting 
we find typical Hi column densities of the BLAs of log 
iV(Hi) « 13 and typical b values of 40 - 100 kms'^ As 
a visual inspection suggests, only 321 of these BLAs are 
single-component, Gaussian shaped absorbers, while the 
remaining 449 systems show evidence for sub-component 
structure and multiple unresolved velocity components. 
Our simulation implies that BLAs predominantly trace 
shock-heated gas in the warm-hot intergalactic medium 



p. Richter, T. Fang, and G.L. Bryan: Simulations of BLAs at low redshift 



11 



0) 



3 

Z 




10 20 30 40 50 
velocity offset [km s~>] 



2 
1.8 



> 1.6 
o 

60 1.4 

o 



1.2 
1 



1 


1 1 1 1 

• 


1 








7- 


. .•%.. •» 


• 


_ # • 






• 

• a 

. *• 
• • 


• '^/^ • 


» - 






















1 1 1 1 


1 


1.6 


1.8 2 2.2 


2.4 




log bni 





Fig. 8. Upper panel: the distribution of velocity offsets 
between broad Lja absorption and associated Ovi ab- 
sorption is shown (gray line = total sample; black line = 
high-quahty sample). Lower panel: O VI 6 values are plot- 
ted against Hi b values (6 in units kms~^). The solid line 
indicates 6(0 vi) = 6(Hi)/4, as expected for pure ther- 
mal broadening (gray dots = total sample; black dots = 
high-quality sample). 



(WHIM) in the temperature range log T = 4.4 — 6.2. 
However, photoionized systems at lower temperatures (log 
T < 4.3) also contribute to the population of BLAs (^ 27 
percent by number) . BLAs that belong to the WHIM typi- 
cally have overdensities of 20 and characteristic ionization 
fractions of log /h ~ 5.6. At these densities, both col- 
lisional ionization and photoionization contribute to the 
total ionization fraction in the gas. We find that in the 
characteristic BLA density range log /h scales linearly 
with log T, so that for pure thermal line broadening the 
ionization fraction can be estimated from the BLA line 
width. The total hydrogen column densities of the BLAs 
typically lie in the range log A^(Hii) = 18 — 20 with a 
peak near log A''(Hii) = 18.8. Most of the photoionized 
systems have logarithmic total hydrogen column densities 
< 18.5 and thus do not contribute significantly to the the 
total baryon content of BLAs. We find that non-thermal 
line broadening mechanisms contribute (on average) with 
~ 10 percent to the total b value of a WHIM BLA, sug- 



gesting that for most systems the measured line width 
represents a good indicator for the thermal state of the 
absorbing gas. Therefore, our simulation suggests that it 
is possible to estimate the total gas column density in 
BLAs from the measured b values and Hi column densi- 
ties and to estimate the baryon content of these systems 
in observational data. Our data imply that this indirect 
method works well for systems with total logarithmic hy- 
drogen column densities < 20.5 and for a BLA sample that 
is sufficiently large. The total baryon content of all BLAs 
in our simulation in the range log A'^(Hll) < 20.5, b > AO 
kms-i and log (N/b) > 10.7 (754 systems; dAf/dz « 38) 
is r2b(BLA) = 0.0121 /i65 which represents ~ 25 per- 
cent of the total baryon budget in our simulation. This 
result supports the idea that BLAs represent a significant 
baryon reservoir at low redshift. If we restrict our total 
BLA sample to the range log 7V(Hii) < 20.5, 6 = 40 - 100 
kms~-^ and log (N/b) > 11.3, we obtain {dAf / dz)BhA ~ 18 
and ^^(BLA) = 0.0058 /les in agreement with recent 
FUV observations. 

Our simulation indicates that future observational 
studies of BLAs at low z will be of great importance to 
study the baryon content of the WHIM, but for a reliable 
estimate of rif,(BLA) a large sample (n > 10) of high- 
resolution, high S/N FUV spectra of low redshift QSOs 
will be required. The Cosmic Origin Spectrograph (COS), 
originally scheduled to be implemented on the Hubble 
Space Telescope (HST) in 2004, would have been the ideal 
instrument for these observations, but in view of the un- 
certain status of the COS mission and other future FUV 
instruments it remains unclear whether such data will be 
available in the near future. Next to observations, addi- 
tional simulations of the low-redshift IGM with higher res- 
olution and more realistic physics (e.g., galaxy feedback, 
metal-line cooling, non-equilibrium processes, etc.) will be 
important to investigate the physical conditions of BLAs 
in more detail and to re-examine the results presented in 
this paper. 

Acknowledgements. P.R. acknowledges financial support by 
the German Deutsche Forschungsgemeinschaft, DFG, through 
Emmy-Noether grant Ri 1124/3-1. T.F. is supported by 
the NASA through Chandra Postdoctoral Fellowship Award 
Number PF3-40030 issued by the Chandra X-ray Observatory 
Center, which is operated by the Smithsonian Astrophysical 
Observatory for and on behalf of the NASA under contract 
NAS 8-39073. G.L.B. is partially supported by NSF grant 
AST0507161. We thank Chris McKee for helpful comments. 

References 

Bryan, G. 1996, Ph.D. thesis, Univ. Illinois, Urbana- 

Champaign 
Gen, R., & Ostriker, J. 1999, ApJ, 514, 1 
Danforth, C.W., & ShuU, J.M. 2005, ApJ, 624, 555 
Dave, R., Hernquist, L., Weinberg, D.H., & Katz, N. 1997 ApJ, 

477, 21 

Dave, R., Gen, R., Ostriker, J., et al.2001, ApJ, 552, 473 
Dave, R., & Tripp, T.M. 2001, ApJ, 553, 528 
Eisenstein, D.J., & Hu, W. 1998, ApJ, 498, 137 



12 



P. Richter, T. Fang, and G.L. Bryan: Simulations of BLAs at low redshift 



Fang, T., & Bryan, G.L. 2001, ApJ, 561, L31 

Fang, T., Marshall, H.L., Lee, J.C., Davis, D.S., & Canizares, 

C.R. 2002, ApJ 572, L127 
Ferland, G.J., Korista, K.T., Vcrncr, D.A., Fergueson, J.W., 

Kingdon, J.B., & Vcrner, E.M. 1998, PASP, 110, 761 
Haardt, F., & Madau, P. 1996, ApJ, 461, 20 
Nicastro, F., et al.2005. Nature, 433, 495 
Norman, M.L, & Bryan, G. 1999, in Numerical Astrophysics, 

ed. S.M. Miyama, K. Tomisaka, & T. Hanawa (ASSL 

Vol.240, Dordrecht: Kluwcr), 19 
Richter, P., Savage, B.D., Tripp, T.M., & Sembach, K.R. 2004, 

ApJS, 153, 165 

Richter, P., Savage, B.D., Sembach, K.R., & Tripp, T.M. 2006, 

A&A, 445, 827 

Savage, B.D., Sembach, K.R., Tripp, T.M., & Richter, P. 2002, 
ApJ, 564, 631 

Sembach, K.R., Tripp, T.M., Savage, B.D., & Richter, P. 2004, 

ApJS, 155, 351 
Sutherland, R.S., & Dopita, M.A. 1993, ApJS, 88, 253 
Tripp, T.M., Savage, B.D., & Jenkins, E.B. 2000, ApJ, 534, LI 
Valageas, P., Schaeffer, R., & Silk J. 2002, A&A, 388, 741 
Williger, G.M., Heap, S.R., Weymann, R.J., et al. 2006, ApJ, 

636, 631 

Zhang, Y., Anninos, P., Norman, M.L., & Meiksin, A. 1997, 
ApJ, 485, 496 



